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We study pairwise thermal entanglement in three-qubit Heisenberg models and obtain analytic 
expressions for the concurrence. We find that thermal entanglement is absent from both the anti- 
ferromagnetic XXZ model, and the ferromagnetic XX Z model with anisotropy parameter A > 1. 
Conditions for the existence of thermal entanglement are discussed in detail, as is the role of degen- 
eracy and the effects of magnetic fields on thermal entanglement and the quantum phase transition. 
Specifically, we find that the magnetic field can induce entanglement in the antiferromagnetic XXX 
model, but cannot induce entanglement in the ferromagnetic XXX model. 



PACS numbers: 03.65.Ud, 03.67.Lx, 75.10.Jm. 
I. INTRODUCTION 

Over the past few years much effort has been put into 
studying the entanglement of multipartite systems both 
qualitatively and quantitatively. Entangled states con- 
stitute a valuable resource in quantum information pro- 
cessing Q. Quite recently, entanglement in quantum 
operations and entanglement in indistinguishable 

fermionic and bosonic systems |5|-|7j have been consid- 
ered. Entanglement in two-qubit states has been well 
studied in the literature, as have various kinds of three- 
qubit entangled states |§| [lC|] . The three-qubit entan- 
gled states have been shown to possess advantages over 
the two-qubit states in quantum teleportation flll| , dense 
coding [[L2| and quantum cloning fl3|| . 

An interesting and natural type of entanglement, ther- 
mal entanglement, was introduced and analysed within 
the Heisenberg XXX @, XX fl§, and XX Z (l6) mod- 
els as well as the Ising model in a magnetic field [ |l7[ . 
The state of the system at thermal equilibrium is rep- 
resented by the density operator p{T) = exp (— ^) /Z, 
where Z =tr[exp (~-^p)] is the partition function, H the 
system Hamiltonian, k is Boltzmann's constant which we 
henceforth take equal to 1, and T the temperature. As 
p(T) represents a thermal state, the entanglement in the 
state is called thermal entanglement jl4j| . A complica- 
tion in the analysis is that, although standard statistical 
physics is characterized by the partition function, deter- 
mined by the eigenvalues of the system, thermal entan- 
glement properties require in addition knowledge of the 
eigenstates. 

The Heisenberg model has been used to simulate a 
quantum computer [fLsf , as well as quantum dots [fLs) , 
nuclear spins [l9| , electronic spins |^] and optical lat- 
tices By suitable coding, the Heisenberg interaction 
alone can be used for quantum computation [^2) . The en- 
tanglement in the ground state of the Heisenberg model 
has been discussed by O'Connor and Wootters [ p3[ . 

In previous studies of thermal entanglement analytical 
results were only available for two-qubit quantum spin 
models. In this paper we analyze the three-qubit case, 



i.e. we consider pairwise thermal entanglement in three- 
qubit Heisenberg models. 

A general 3-qubit Heisenberg XYZ model in a non- 
uniform magnetic field B is given by: 

H = HxYZ + -ffmag 

n ( Jl ^^ _i_ J2 „y„v i J* z z \ 

n=l V - / 

3 

#mag = ^ B n q Z n . (1) 
n=l 

We use the standard notation, detailed later, and assume 
a periodic boundary, identifying the subscript 4 with 1 
in the above expressions. For the 3-qubit case even this 
most general scenario is susceptible to numerical analy- 
sis. However, in this paper we shall restrict ourselves to 
special cases of Eq. (Q) for which we are able to provide a 
succinct analytic treatment. 

The 3-site Heisenberg models we will study in this pa- 
per are the following: 

1. The XX model, corresponding to J\ = J 2 , J3 = 
and B = 0. 

2. The XX Z model, for which J x = J 2 , J 3 ^ and 
B = 0. 

3. The XXZ model with uniform magnetic field 
(B\ = B 2 = B 3 ). 

We start in Sec. II by examining the three-qubit XX 
model. In Sec. Ill, IV, and V, we study thermal en- 
tanglement in the XX model, the XXZ model and the 
XXZ model in a magnetic field, respectively. 

During the course of the analysis it will become clear 
that degeneracy plays an important role in thermal en- 
tanglement, as does the presence of magnetic fields. We 
find the critical temperatures involved in the quantum 
phase transition associated with the existence of entan- 
glement in these quantum spin models. 
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II. THREE-QUBIT XX MODEL AND ITS 
SOLUTION 

The three-qubit XX model is described by the Hamil- 
tonian Q 

J 3 

Hxx = 2" 1] + 

n=l 
3 

= J H + a n a n+l) ' ( 2 ) 

n=l 

where a™ (a — x, y, z) are the Pauli matrices of the rt-th 
qubit, <r„ = I (cr^ ± io*) the raising and lowering oper- 
ators, and J is the exchange interaction constant. Pos- 
itive (negative) J corresponds to the antifcrromagnetic 
(ferromagnetic) case. As signalled above, we adopt pe- 
riodic boundary conditions; erf = of , erf = of. We are 
therefore considering a three-qubit Heisenberg ring. The 
XX model was intensively investigated in 1960 by Lieb, 
Schultz, and Mattis fgj]. More recently the XX model 
has been realized in the quantum-Hall system |25[, the 
cavity QED system [^6| and quantum dot spins |[27| for 
a quantum computer. 

In order to study thermal entanglement, the first step 
is to obtain all the eigenvalues and eigenstates of the 
Hamiltonian Eq.(|^). The eigenvalues themselves do not 
suffice to calculate the entanglement. The eigenvalue 
problem of the XX model can be exactly solved by the 
Jordan- Wigner transformation (28). In the three-qubit 
case the eigenvalues are more simply obtained as p"5[ 

E = E 7 = 0, 

Ei = E% = E4 = E5 = — J, 

E 3 = E 6 = 2 J. (3) 

and the corresponding eigenstates are explicitly given by 



IV-o) 


= 1000), 


1^1) 


= 3~ 1/2 (g|001) + <7 2 |010) + |100» , 


|V>2) 


= 3- 1/2 (<? 2 |001) + g|010) + |100» , 


|V>3> 


= 3~ 1/2 (|001) + |010) + |100)), 


m 


= Z- 1 ' 2 (g|110) +(Z 2 |101) + |011)), 




= 3~ 1/2 (g 2 |110) +g|101) + |011» , 


l</<6> 


= 3" 1/2 (|110) + |101) + |011)), 




= |111). 



with q = exp (i27r/3) satisfying 

c/ = l, 

q 2 + q + 1 = 0. (5) 

This set (^) of three-qubit states is itself interesting. Ra- 
jagopal and Rendell jl0| have considered a similar set of 



three-qubit states which they have classified by means of 
permutation symmetries. Here the states |^o), 1^3) > \4>6)i 
and \ip-r) are symmetric in the permutation of any pair 
of particles. We define a cyclic shift operator P by its 
action on the basis \ijk) p9[ 

P\ijk) = \kij). (6) 

Obviously the four states |^o)> | "03 ) i l^te), and I ^7) are 
the eigenstates of P with eigenvalue 1. The other four 
states in the set ^ are also eigenstates of P as follows: 

P\ipi) = q 2 \i>i) (< = 1,4), (7a) 
P\i> j ) = q\*p j ) (3=2,5). (7b) 

This is not surprising since the Hamitonian Hxx as well 
as the other Hamiltonians considered later are invariant 
under the cyclic shift operator. 

For J > (J < 0) the ground state is four (two)- 
fold degenerate. We will see that the degeneracy of the 
system influences thermal entanglement greatly. There 
is no pairwise entanglement in the eigenstate |^o) an d 
Pairwise entanglement exists in the state (i — 
1, 2, 6) and the concurrence between any two different 
qubits is given by 2/3 §@. 

III. THERMAL ENTANGLEMENT IN THE XX 
MODEL 

We first recall the definition of concurrence ||l|] be- 
tween a pair of qubits. Let P12 be the density matrix of 
the pair and it can be either pure or mixed. The con- 
currence corresponding to the density matrix is defined 
as 

C12 - max {Ai - A 2 - A 3 - A 4 , 0} , (8) 

where the quantities A^ are the square roots of the eigen- 
values of the operator 

012 = Pii{oi <8> erf )p* 12 ( (T l ^I) (9) 

in descending order. The eigenvalues of Q12 are real and 
non-negative even though Q12 is not necessarily Hermi- 
tian. The values of the concurrence range from zero, for 
an unentangled state, to one, for a maximally entangled 
state. 

The state at thermal equilibrium is described by the 
density matrix 

p(T) = I exp (-0H) , 
1 7 

= -;^exp(-^ fe )|^)(Vfe|- (10) 

fc=0 

where j3 = 1/T. From Eq.(^|), the partition function is 
obtained as 
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2 + 4e /3J + 2e~ 2/3J . 



(11) 



/From Eqs.([3]) and (|l0|), the density matrix can be writ- 
ten as 



P(T) = |[K>o><^> 



(^7 1 

e^d^iXV'il + l^)^! 

I^K^I + I^K^I) 



+e" 



-2/3J 



(1^3)^31 + 1^6)^61)]. 



(12) 



In this paper we consider only pairwisc thermal entan- 
glement, and so we need to calculate the reduced density 
matrix p\2{T) = tr 3 (p(T)). We denote the reduced den- 
sity matrix tT 3 [\ip il )(ip il \ + ... + \ip iN )(iJj iN \] by p^ 2 '"'' x 
From Eq.(0), we obtain 
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(13a) 



(13b) 



(13c) 



(13d) 



(13e) 



The last two reduced density matrices will be used later. 
/From Eqs. (|l2|) and ( |l^|l3c] ), we obtain 



n m- 1 L (07) + P P J n {12i ^ 4_ p -2/3J„(36)\ 
Pl2(,J j - \P\2 + e Pl2 + e Pl2 ) 



2 

3Z 



with 



/ v 

w y 

y w 

\ v 



3 3/ 1 
- + eJ iJ + -e 



(14) 



2/3J 



2 

2e? J 



+ e" 



-2/3J 



2 

-20J 

i 



(15) 



The square roots of the four eigenvalues of the operator 
£>12 are 



Ai 



2(w-y) _ 2(w + y) 

-, *2 — 



3Z 



3Z 



x -x - 2v 

A 3 - A 4 - — . 



(16) 



iFrom EqsM, (0), and we obtain the con- 

currence \U3\ 



C = — max(|y| - v,0), 



max 



2|e" 



3-2e x 



3(1 + 26^+6- 



-,0 



(17) 
(18) 



where x = (3 J — J/T. The concurrence depends only 
on the ratio of J and T. Due to symmetry under cyclic 
shifts, the value of the concurrence does not depend on 
the choice of the pair of qubits. 

^From we see that entanglement appears only 
when 



2\z~ 2 - z\ - 3 - 2z - z~ 



3(1 

or in other words 



2z- 



>0, 



2\z~ 2 - z\ - 3 - 2z - z~ 2 > 0, 



(19) 



(20) 



where z = exp(x). We now consider two different cases: 
Case 1. Antiferromagnetic system; J > 0;z~ 2 — z < 0. 
In this case relation (EG) requires 



z- 2 < 



-1 



(21) 



which is impossible. So there is no entanglement when 
J > 0. 

Case 2. Ferromagnetic system; J < 0,z~ 2 — z > 0. 
Relation (GOT) becomes 



or 



4z - 3 > 0. 



f{z) = Az 3 + 3z 2 - 1< 0. 



(22) 



(23) 



The function f(z) is an increasing function of the positive 
real argument z and relation (|23|) is valid iff < z < z c , 
where the critical value z c determined by f(z c ) = is 
0.4554; that is, x < -0.7866. For fixed J, we obtain the 
critical temperature T c — 1.21736| J|, above which there 
is no thermal entanglement. The critical temperature 
depends linearly on the absolute value of J. 
In the ferromagnetic case the concurrence 



1 - 4z 3 - 3z 2 
3(1 + 2z 3 + z 2 ) 



,0 



(24) 



reaches its maximum value of 1/3 when z — > 0, that is 
when x — > — oo. Since the entanglement is a monotonic 
increasing function of C this means that the entanglement 
attains its maximum value for zero temperature, when 



3 



J is finite and nonzero. For finite temperatures, this 
maximum is also attained when J — > — oo. 
In summary, we find that 

Theorem 1: The XX model is thermally entangled if 
and only if J < — 0.7866T; maximum entanglement is 
attained when T — ► or J — > — oo. 

The above discussion shows that in our 3-qubit model 
pairwise thermal entanglement occurs only in the ferro- 
magnetic case. This result differs from that for the two- 
qubit XX model, for which thermal entanglement exists 
in both the antiferromagnetic and ferromagnetic cases 

& 

For the ferromagnetic case the states |i/ , 3) and It/jr) con- 
stitute a doubly-degenerate ground state. Eq.(g4|) shows 
that the concurrence C = 1/3 at zero temperature. As 
noted in the last section the concurrence for any two 
qubits in the state ^3) or \tpo) is 2/3. Here the value 1/3 
appears due to the degeneracy. In fact, at zero temper- 
ature, the thermal entanglement can be calculated from 

Pia^COi)' After normalization it is easy to obtain the 
concurrence, which is just 1/3. 



IV. THE ANISOTROPIC HEISENBERG XXZ 
MODEL 

We now consider a more general Heisenberg model, the 
anisotropic Heisenberg XXZ model, which is described 
by the Hamiltonian |32| 



H 



AJ 



xxz - H X x + —z- ^{KK+i - 1) = 



(25) 



As in the last section, since Z is always positive, we 
need only consider 



/(A,z) = |y|-« 



z 2A \z- 2 -z\---z 2A+1 --z 2A - 2 
1 1 2 2 



(28) 



to determine whether entanglement occurs or not. Again, 
we have to consider two different cases: 

Case 1. When J > (z > 1), namely the antifer- 
romagnetic XXZ model, the condition on /(A, z) leads 
to 



z 



2A-2 



2^(A-1) 



< -1, 



(29) 



which is impossible. So there is no entanglement in this 
case, irrespective of A. 

Case 2. When J < (z < 1), namely the ferromag- 
netic XXZ model, the condition /(A, z) > gives 



z 2A-2 _ 4z 2A+l _ 3 > 0. 



(30) 



We consider some special values of A. 

(1) A > 1: For A = 1 the relation @ implies z 3 < 
— 1/2 which is impossible. So there is no entanglement 
in the XXX model. We can further prove that there is 
no entanglement for A > 1. In fact, it is easy to see that 



/(A,z) < z 2[A -^ -3 < 0, 



(31) 



where we have used the inequalities z 2A+1 > and 
z 2(A-i) < 1 f or A > 1 and z < 1. This means C = 
and thus there is no entanglement. 

(2) A = 1/2: In this case 33 1 the entanglement condi- 
tion is obtained as 



where A is the anisotropy parameter. The model reduces 
to the XX model when A = 0, and the isotropic Heisen- 
berg XXX model when A = 1. 

It is straightforward to check that the added 
anisotropic term Hxxz ~ Hxx commutes with Hxx- 
Therefore the eigenstates of the XXZ model are still 
given by Eq. ([l]), now with the different eigenvalues 

E Q = E 7 = 0, 

E x = E 2 = E 4 = E 5 = -2J(A + i), 



4z d 



3z - 1< 0, 



£ 3 = E e = -2 J(A - 1) 



(26) 



Following the procedure of the previous section, we ob- 
tain the concurrence, which is of the same form as Eq.jl^) 
with however the parameters v, w, y, and the partition 
function Z now given by 



3 1 

v = 7; + t: z 



2 A 



(22- 

-2\ 



2 2 

w = z 2A (2z + z 

2Ar -2 \ 

y = z (z -z), 

Z = 2 + 2z 2A (2z + z~ 2 ). 



2 ). 



(27) 



(32) 



which is an increasing function of z. So the model is 
entangled iff < z < z c w 0.298, where z c is determined 
as a root of 4z 3 + 3z — 1 = 0. 

(3) A = —1/2: This is an interesting case whose 
importance has been emphasized recently [p^| . ^From 
the eigenvalues we see that the excited state of the sys- 
tem is 6-fold degenerate when A = —1/2. The func- 
tion /(A, z) now reduces to z~ 3 — 7, from which the 
critical values z c and T c are obtained analytically as 
z c = 7~ 1 / 3 ,T C = 3|J|/ln7w 1.5417|J|. 

(4) The limit case A — > — 00: The critical value z c is 
now determined by z~ 2 — 4z = 0, i.e., z c = 4 -1 / 3 . There- 
fore the critical temperature T c — 3| J|/ln4 rj 2.164|J|. 

Finally, for more general values of the anisotropy pa- 
rameter we need to resort to numerical calculations. 
Fig.l is a plot of the critical temperature as a function of 
the anisotropy parameter A. From this we see that the 
critical temperature decreases as A increases, and reaches 
the asymptotic value T c = 2.1640| J\ as A — > —00. 
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as indicated in Fig. 2. 




f(A) 



-7 -6 -5 -4 



FIG. 2. The function /(A, z) with respect A for 
z = 0.6295, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1 (from left to right). 



FIG. 1. The critical temperature T c as a function of the 
anisotropy parameter A. The exchange constant J = — 1. 

We now give further analytical results for the case A < 
1 and z < 1. Consider /(A, z) as a function of A. Then, 
from 



df(A,z) 
dA 



= (In z)z^(z- 2 -4z) 

= 0, when z = zq = 4~ 1//3 
> 0, when z > zq; 
< 0, when z < Zq. 



0.62996; 



(33) 



we see that /(A, z) is an increasing (decreasing) function 
when z > zq (z < zq). We consider these cases separately. 

Case 2a. When z = zq, /(A, z ) = — 3 < 0. So there 
is no entanglement in this case. 

Case 2b. When z > Zq, the function /(A,z) is an 
increasing function which reaches its maximum when 
A — > 1. Since we have seen that there is no entangle- 
ment when A = 1 



f(A,z) < f{l,z) < forz>z , 



(34) 



which means that there is no entanglement when z > zq. 

Case 2c. The case z < zq. Define the z-dcpendent 
point A z by f(A z , z) = where 



A, 



2(3J 



Tn 



-Az 



< 1. 



(35) 



Thus from Eq.(p3|) we know that f(A,z) > when 
A < A z for all z < zq, which is just the condition for 
entanglement. 

In Fig. 2 we give plots of /(A, z) for z = 
0.6295,0.6,0.5,0.4,0.3,0.2,0.1. Note that A z is a de- 
creasing function of z and that 



A, 



— oo when z 
1 when z 



z 0'i 

0, 



(36) 



In summary, we have 

Theorem 2 The XXZ model exhibits thermal entan- 
glement only when z < z$, that is, J < —.462 IT, and 
A < A z . 

Note that Theorem 2 is entirely consistent with The- 
orem 1, since A z = in Eq.(35) corresponds to J = 
-0.7866T. 

In Fig. 3 we plot the concurrence as a function of the 
anisotropy parameter A and exchange constant J. The 
figure shows that there is no thermal entanglement for 
the antiferromagnetic (J > 0) XXZ model, nor for the 
ferromagnetic (J < 0) XXZ model when A > 1. 
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FIG. 3. The concurrence as a function of A and J. The 
temperature T = 1. 

To end this section we investigate the concurrence at 
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zero temperature. Nonzero concurrence occurs for the 
case A < 1, J < 0. In this case, the doubly-degenerate 
ground state consists of ^3) and j?/^). We may calculate 
the concurrence C = 1/3 directly from the density matrix 

J36) 
Pl2 ■ 



V. EFFECTS OF MAGNETIC FIELDS 

In this section we consider the effect of magnetic fields 
on thermal entanglement. The XX Z model with uniform 
magnetic field B along the z direction is given by 



H 



xxzm — Hxxz 



(37) 



It is easy to check that the added magnetic term com- 
mutes with the Hamiltonian Hxxz- Therefore the eigen- 
states of the XX Z model are given by Eq.(Q). The eigen- 
values are now 

E Q = -W 

E 1 =E 2 = -2J(A + ^)-B 

E 3 = -2J(A-l)-B, 

E 4 = E 5 = -2J(A + ^) + B, 

E 6 = -2J(A-1) + B. 

E 7 = 3B (38) 

We see that the magnetic field partly removes the degen- 
eracy. 

With a derivation completely analogous to that of Sec. 
Ill and Sec. IV, the reduced density operator is 



P12 



u 

_2_ [ w y 

3Z I y w 

v 



(39) 



with 



v = ^e-W B + ±e-P B z 2A (2z + z~ 2 ) , 

w = cosh(/3 J B)z 2A (2z + z~ 2 ) , 
y = cosh((3B)z 2A (z- 2 - z) , 
Z = 2cosh(3/3 J B) + 2 cosh(/3 J B)z 2A 
x(2z + z~ 2 ) 

The concurrence is then given by 
4 

C = — max(|y| - ^/uv, 0). 



(40) 



(41) 



As an immediate consequence we see that the concur- 
rence is an even function of the magnetic field. 



As the quantities Z, u, v are all positive, for con- 
venience we consider the quantity y 2 — uv instead of 
\y\ — y/uv. Thermal entanglement occurs when 



V 



uv 



h cosh(2/3£) -g > 0, 



(42) 



where 



.9= J[9 + ^ 4(A - 1) (2z 6 + 8^-l)], 

h = I Z 2A [ z 2A (z -2 _ z) 2 _ ( g z + 3z _aj] _ (43) 

We now consider the effect of a magnetic field on the 
thermal entanglement. 

We first consider the XXX model, A = 1, which does 
not exhibit thermal entanglement when B = 0. One 
might expect that the magnetic field would induce ther- 
mal entanglement. It is easy to see that 

2{y 2 - uv) = cosh(2/3B) (z 6 - 8z 3 - 2) - (z 3 + 2) 2 . 

(44) 

If z < (4 + 3\/2) 1/3 ~ 2.02, z 6 - 8z 3 - 2 < and thus 
y 2 — uv < for any B. So for this range of z values 
there is no thermal entanglement no matter how strong 
the magnetic field is. However, when z > (4 + 3^2) 1 ^ 3 , 
z 6 — 8z 3 — 2 > and the condition for entanglement 
becomes 



cosh(2/3B) > 



(z 3 + 2) 2 
6 - 8z 3 - 



(45) 



which can be fulfilled for strong enough B. So a magnetic 
field can induce entanglement in the XXX model when 
z > (4 + 3\/2) 1 / 3 . 

Now consider the case A = —1/2. ^From Eq.(|43|) we 
obtain 



h =\(P 2 -5p-5), 
.9= J(ll + 8p-p 2 ), 
h-9= i(3p 2 -18p-21), 



(46) 
(47) 
(48) 



which are parabolas in p 



5 , as shown in Fig. 4. We 



consider three different cases: 

Case 1: p < pi = 5/2 + 3\/5/2, In this case h < 
0,g > 0, h — g < and y 2 — uv < 0. So there is no 
thermal entanglement. 

Case 2: p\ < p < p 2 = 7. In this case h > 0, g > 
0, h — g < 0. So y 2 ~ uv > 0, and so entanglement 
appears if the magnetic field is strong enough. 

Case 3: p 2 < P- In this case h > 0, h—g > and y 2 —uv 
is always positive; that is, here the XXZ model exhibits 
thermal entanglement for any magnetic field. Note that 
p 2 = z~ 3 where z c is the critical value given in last sec- 
tion. 
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The above two models show that the magnetic field can 
either induce entanglement in a non-entangled system or 
extend the entanglement range for an already entangled 
system. 



FIG. 5. Concurrence as a function of the magnetic field 
B and the exchange constant J. The temperature T = 1 and 
the anisotropy parameter A = 1. 
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FIG. 4. The functions h, g, and h — g in terms of p = z . 

In Fig. 5 we plot the concurrence as a function of the 
magnetic field B and exchange constant J. At B = 
there is no thermal entanglement. The entanglement in- 
creases with the magnetic field \B\ until it reaches a maxi- 
mum value, then decreases and gradually disappears. We 
can clearly see that there is no thermal entanglement for 
the ferromagnetic case, while thermal entanglement ex- 
ists for the antiferromagnetic case. In other words, we 
can induce entanglement in the antiferromagnetic XXX 
system by introducing a magnetic field, but cannot in- 
duce entanglement in the ferromagnetic XXX system 
for any strength of magnetic field. 
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FIG. 6. Concurrence as a function of T for different mag- 
netic fields B — l(solid line), 3/2(dashed line), and 2(circle 
point line). 

Fig. 6 gives a plot of the concurrence as a function of 
the temperature for different magnetic fields. One can 
see that there exist critical temperatures above which 
the entanglement vanishes. It is also noteworthy that 
the critical temperature increases as the magnetic field 
B increases. Consider the interesting case B = 2. We 
observe that the concurrence is zero at zero temperature 
and there is a maximum value of concurrence at a finite 
temperature. The entanglement can be increased by in- 
creasing the temperature. The maximum value is due 
to the optimal mixing of all eigcnstates in the system. 
When considering zero temperature we find that there 
arc different limits for different magnetic fields. Actually 
a more general result exists 

lhnC(A,5,l,T) = i for A > \B\ - 1/2. 

= ^ for A = \B\ - 1/2, 

= for A < \B\ - 1/2. 

The special point T = 0, A = B - 1/2 (B > 
is assumed without loss of generality), at which the en- 
tanglement undergoes a sudden change with adjustment 
of the parameters A and B, is the point of quantum 
phase transition ]3q| . The quantum phase transition 
takes place at zero temperature due to the variation of 
interaction terms in the Hamiltonian. By examining the 
eigenvalues (|3^) we can understand the phase transi- 
tion. When A = B — 1/2, the ground state contains 
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the three-fold degenerate states \ipo), and \tp2)- One 
may calculate the thermal entanglement from the den- 



sity matrix p^ 12 '*(13d) and find the concurrence to be 



hospitality of the Laboratoire de Physique Theorique des 
Liquides, Paris University VI. 



2/9. When A > B — 1/2, the ground state contains 
the two-fold degenerate states and |-02)- The con- 
currence has the value 1/3 as calculated from p^i^c). 



When A < B — 1/2, the ground state is li/'o) and not 
degenerate. And the concurrence is zero in this case. 



VI. CONCLUSIONS 

Apart from being a fundamental property of quantum 
mechanics, it appears that entanglement may provide an 
important resource in quantum information processes. 
One source of entanglement is provided by magnetic sys- 
tems, such as those modelled in this paper. Within the 
current state of knowledge, only measures for pairwise 
entanglement are available. Thus, in order to study the 
entanglement properties of systems more complex than 
those simply involving two qubits, it is necessary to adopt 
a procedure whereby one traces out a subsystem, leav- 
ing effectively only a two-qubit system for which we can 
calculate the concurrence, which in turn gives a mea- 
sure of the entanglement. Using this procedure, we have 
studied pairwise thermal entanglement in the following 
Heisenberg models; the XX model, the XX Z model and 
the XXZ model in a magnetic field. We obtained ana- 
lytical expressions for the concurrence, which indicated 
no thermal entanglement for the antiferromagnetic XXZ 
model, nor for the ferromagnetic XXZ model when the 
anisotropy parameter A > 1. Conditions for the exis- 
tence of thermal entanglement were studied in detail. 
The effects of magnetic fields on entanglement were also 
considered. We found that the magnetic field can induce 
entanglement in the antiferromagnetic XXX model, but 
cannot induce entanglement in the ferromagnetic XXX 
model, no matter how strong the magnetic field is. 

In this paper we have extended previous work on ther- 
mal entanglement from two-qubit models to three qubit 
models, concentrating on those systems where the pair- 
wise entanglement can be studied analytically. It would 
be an attractive proposition to extend further the inves- 
tigation of such Heisenberg models to the A-qubit case, 
which are under consideration. 
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